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Abstract 

Charmonium production at heavy-ion colliders is considered within the comover interaction 
model. The formalism includes both comover dissociation of J/V''s and possible secondary J/^ 
production through recombination. The estimation of this effect is made without involving free 
parameters. The comover interaction model also incorporates a comprehensive treatment of initial- 
state nuclear shadowing. With these tools, the model properly describes the centrality, transverse 
momentum and rapidity dependence of experimental data from PbPb collisions at the LHC energy 
of ^/s = 2.76 TeV. We present predictions for PbPb collisions at ^/s = 5.5 TeV. 
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1 Introduction 



Lattice QCD calculations predict that, at sufficiently large energy densities, hadronic matter undergoes 
a phase transition to a plasma of deconfined quarks and gluons (QGP). Substantial activity has been 
dedicated to the research of high-energy heavy-ion collisions in order to reveal the existence of this 
phase transition and to analyze the properties of strongly interacting matter in the new phase. The 
study of quarkonium production and suppression is among the most interesting investigations in this 
field since, in the presence of a QGP, the charmonium yield would be further suppressed due to color 
Debye screening [1]. Indeed, such an anomalous suppression was first observed in PbPb collisions at 
top CERN SPS energy [2]. Alternatively, the SPS experimental results could also be described in 
terms of final state interactions of the cc pairs with the dense medium created in the collision, the so- 
called comover interaction model (CIM) [3llH[5l|6]. This model does not assume thermal equilibrium 
and, thus, does not use thermodynamical concepts. 

The theoretical extrapolations to RHIC and LHC energies were guided mainly by two trends. 
On the one hand, models assuming a deconfined phase during the collision pointed to the growing 
importance of secondary J /ijj production due to regeneration of cc pairs in the plasma. The total 
amount of cc pairs is assumed to be created in hard interactions during the early stages of the collision. 
Then, either using kinetic theory and solving rate equations for the subsequent dissociation and 
recombination of charmonium [7l[8l|9], or assuming statistical coalescence at freeze-out |10 1 H H [T2]. 
one obtains the final J /tp yield. On the other hand, the CIM with only dissociation of J/ip^s predicted 
[13] a stronger suppression at RHIC than at SPS due to a larger density of produced soft particles 
in the collision. It also predicted a stronger suppression at mid-rapidity -where the comover density 
is maximal" than at forward rapidities. Nevertheless, measurements of J/ip production in AuAu 
collisions at RHIC, y/s = 200 GeV, gave surprising results: the suppression at mid-rapidity was on 
the same level as at SPS I15|. Furthermore, the suppression at forward rapidity in AuAu collisions 
was stronger than at mid-rapidity. 

These facts showed that in the CIM, which is based on the well-known gain and loss differential 
equations in transport theory, the introduction of a recombination term is actually required for a 
detailed balance |16| . In the present work we use the updated version of the CIM |T6j, that allows 
recombination of cc pairs into secondary J/tp^s. We will estimate this effect using the density of 
charm in proton-proton collisions at the same energy. Therefore, the model does not involve any 
additional parameters. In Section [2] we present the details of the model; the effects related to the 
initial-state shadowing and nuclear absorption, together with comover dissociation and regeneration 
are described. In Section|3]we present our results for PbPb collisions at the LHC energy of = 2.76 
TeV and predictions at ^/s = 5.5 TeV. Conclusions and final remarks are given in Section [4] 
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2 Description of the model 



In this section, we give a short and updated description of the CIM. The present version of this 
model contains a comprehensive treatment of initial-state nuclear effects, such as nuclear shadowing, 
together with nuclear absorption, interaction with the co-moving matter and the recombination of cc 
into secondary J/ip^s. 

Nuclear effects in nucleus-nucleus collisions are usually expressed through the so-called nuclear 
modification factor, R'^^{b), defined as the ratio of the J ftp yield in AA and pp scaled by the number 
of binary nucleon-nucleon collisions, n{b). We have then 



n(6)div//^/dy 
/ dh aMb) nib, s) S-'^b, s) Sf^^{b, s) S^"{b, s) 



1) 



/ d'^saAB{b)n{b,s) 

where (TAB{b) = 1 — exp[—appABTAB{b)], TAsib) = J d'^sTA{s)TB{b — s) is the nuclear overlapping 
function and T^(6) is obtained from Woods-Saxon nuclear densities [17J. In eq. ([T]), 

n{b,s) = appABTA{s)TB{b- s)/aAB{b) , (2) 

where upon integration over d^s we obtain the number of binary nucleon-nucleon collisions at impact 
parameter b, n{b). 

The three additional factors in the numerator of eq. ([T]), S"^'*, S^^ and 5^^°, denote the effects 
of nuclear absorption, shadowing and interaction with the co-moving matter -both dissociation and 
recombination-, respectively. They will be defined below. 

2.1 Nuclear absorption 

At low energies the primordial spectrum of particles created in scattering off a nucleus is mainly 
altered by interactions with the nuclear matter -spectators nucleons- they traverse on the way out to 
the detector. This effect is called nuclear absorption and is usually parameterized within a probabilistic 
Glauber model. One has 

gabs/, X _ [1 - exp{-A Ta{s) aabsW - exp{-B Tsjb - s)aabs)] , . 

al,,ABTA(s)Ts{b-s) ' 

Several novel features are expected to appear at high energies due to the change in the space-time 
picture of particle production and the coherence effects [18]. In this way, the break-up cross section 
becomes aahs = at high energies |19^ [20l I21j . From now on, we will take S""^^ = 1 for LHC energies, 
which implies a negligible nuclear absorption. 
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2.2 Shadowing 

Coherence effects wiil iead to nuclear sliadowing for both soft and hard processes at high energy, and 
therefore also for the production of heavy flavor. Shadowing can be calculated within the Glauber- 
Gribov theory [22J making use of the generalized Schwimmer model of multiple scattering [23j . In this 
case the second suppression factor in eq. ([T]) is given by 

^''^^'"'^^ = 1 + AF{yA)TA{s) 1 + BF{yB)TB{h-s) ' 

where the function F{y) encodes the dynamics of shadowing. Following the spirit of the model 
presented in |24l [25] , where shadowing corrections are given without free parameters in terms of the 
triple-Pomeron coupling determined from diffraction data, one can write: 

fVmax \ dCT^^^ 



rVmax 

F{y) = 47r / dy 



ap{y) dydt 



max 

) — exp(Aymj„)] /A , (5) 

t=0 



It represents the coherence effects due to shadowing corrections expressed as the ratio of the triple- 
Pomeron cross section over the single-Pomeron exchange. We take ymin = £n{RA'mN /\/3), where m^v 
is the nucleon mass and Ra = 0.82 + o.58 fm - the gaussian nuclear radius. The value of ymax 
depends on the rapidity of the produced particle, ymax = \'a.{s/mj)/2 it y with the +(— ) sign if the 
particle is produced in the hemisphere of nucleus B[A). niT is the transverse mass of the produced 
particle. We have used A = 0.13 and C = 0.04 fm-^ (C'/A = 0.31 fm^) corresponding to the Pomeron 
intercept ap(0) = 1.13. With the shadowing resulting from the above equations a good description 
of the centrality dependence of charged multiplicities is obtained both at RHIC j[26j and LHC |i27j 
energies. 

The energy dependence of these different effects can be summarized as follows. Particle production 
at SPS is dominated by low-energy effects, i.e. nuclear absorption and small nuclear shadowing, while 
RHIC already belongs to the high-energy regime, where nuclear shadowing becomes relevant, and the 
combined effect of shadowing and energy-momentum conservation should be accounted for at forward 
rapidities. At LHC, shadowing is expected to be strong while nuclear absorption is a small effect and 
can be neglected. We will now proceed with the discussion of comover interaction effects. 

2.3 Dissociation by comover interaction and recombination 

The CIM was developed to explain both the suppression of charmonium yields [HI HI El El [l3l [281 
[29] and the strangeness enhancement |30[ [3T] in nucleus-nucleus collisions at the SPS. Neglecting 
possible recombination effects, the rate equation governing the density of charmonium in the final 
state, Nj/^, can be written in a simple form assuming a pure longitudinal expansion of the system 
and boost invariance. For an AA collision the density of J/ijj at a given transverse coordinate, s, 
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impact parameter b, and rapidity y is given by 



dA^ 



J/V 



dr 



{b,s,y) = -acoN''°{b,s,y)Nj/^,{b,s,y) 



(6) 



where aco is the cross section of charmonium dissociation due to interactions with the co-moving 
medium, with density N'^°. It is found from fits to low-energy experimental data to be aco = 0.65 
mb [5]. To incorporate the effects of recombination, we have to include an additional gain term 
proportional to the squared density of open charm produced in the collision. Then eq. ^ is generalized 
to 



dN 



dr 



ib,s,y) 



-(Tr- 



N"'{b,s,y)Nj/4b,s,y) - N^ib, s,y)N,{b, s,y) 



(7) 



where we have assumed that the effective recombination cross section is equal to the dissociation cross 
sectiorQ This extension of the model therefore does not involve additional parameters. 
Equation Q cannot be solved analytically. We approximate its solution as 



S'°ib,s,y) 



exp 



-(Jr. 



N'°ib,s,y) 



Ncib,s,y)N-c{b,s,y) 



In 



N'°{b,s,y) 



(8) 



Nj/^i,{b,s,y) 

resembling the exact solution of eq. ([g]), since the first term in the exponent of eq. ([s]) is exactly the 
survival probability of a J/^ interacting with comovers |13| . The density of comovers is calculated 
following the same lines as in [27J together with the proper shadowing correction. Then 



N'^%b,s,y) = N'^sib,s,y)S!tib,s,y), 



(9) 



where S^j^ denotes the shadowing for light particles. The non-shadowed (NS) multiplicity of comovers 
is taken as proportional to the number of nucleon-nucleon collision according to eq. ^ 



N'^sib,s,y) = NPP{b,s,y)nib,s), 



(10) 



where N^^ represents the comover density in pp collisions, essentially NPP{y) = KdN"'' / dyY^ . In 
fact, when using at mid-rapidity the value of (dA^^''/dy)^^Q = 3.9, compatible with the inelastic 
value of dN^P /dr] at y/s = 2.76 TeV, a good agreement with experimental data on charged particle 
multiplicities is obtained [27] . 

The density of open and hidden charm in AA collisions, Nc,Nc and Nj/^, respectively, can be 
computed from their densities in pp collisions as N^^{b,s) = n{b, s)S1^Q{b, s)N'§p , with similar ex- 
pression for N-^"^ and -^jy^- Here n{b, s) is given by eq. Q and Sf^q is the shadowing factor for heavy 
quark production, given by eq. Q. Then eq. ([s]) becomes 



S'"'{b,s,y) 



exp 



-(Jr. 



N'°{b,s,y)-C{y)n{b,s)Sfj%{b, 



In 



(11) 



* These two cross-sections have to be similar but not necessarily equal. We have taken the simplest possibility. 
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where 



C{y) = ^ J. ' = ^ J ' . (12) 
dNpy/dy appdapy/dy 



The quantities in the rightmost term in eq. (12) are all related to pp collisions at the corresponding 
energy. The value for ddp/^/dy can be taken from ALICE [3^ or a model for extrapolation of the 
experimental results [33]. The cc pairs are mostly in charmed mesons, such as D and D* and the 
corresponding values could be extracted from the experiment |3l], leading to an estimation of dapp/dy, 
as we will discuss below. For cjpp we use the non-diffractive value app=5A mb [27j. 

With (Tco fixed from experiments at low energy, where recombination effects are negligible, the 
model, formulated above, should be self-consistent at high energie^ 



3 Results 

Our results for the centrality dependence of the J/ip nuclear modification factor in PbPb collisions at 
2.76 TeV are presented in Fig. [T] compared to ALICE experimental data j35ll36t [37 1 l38j at mid and for- 
ward rapidities. The different contributions to J/tp suppression are shown. Note that the initial-state 
effect is just the shadowing, which can induce a suppression of Raa = 0.6 for the more central colli- 
sions [39] . The combined effect of shadowing and comover dissociation gives a too strong suppression 
compared to experimental data. We therefore proceed to estimate the effect of recombination. 

For the charmonium cross section pp measurements around mid-rapidity are available from ALICE 
|32j at 2.76 TeV which corresponds to = 3.73 /xb. We consider that realistic values of C at mid- 

rapidity at 2.76 TeV are in the range of a minimum value of 2 up to a maximal value of 3 which 
corresponds to a cross section ~ 0.6 0.8 mb. This agrees with the estimated values in |40j . and 
corresponds to a cr*g* around 5 mb, which agrees well with experimental data [33]. These values are 
higher than the ones reported in [51], where data is also reproduced. Note nevertheless that there 
is no contradiction, since in ^1] the initial-state shadowing is not considered. This shadowing, that 
affects also heavy fiavours, would imply an extra suppression leading naturally to the choice of higher 
input charm cross sections dapp/dy. In other words, one can consider that in the approach developed 
in [JT] the choice of smaller ~ 0.3 0.4 mb, takes into account an effect of shadowing that 

reduces the input charm cross section up to 1/2. 

We expect the effect of recombination to be stronger at mid than at forward rapidities. At y 7^ 
the recombination term is smaller since the rapidity distribution of D, D* is narrower than the one of 

^Note, however, that CTco could change when the energy increases. We do not expect this effect to be important and, 
since we are unable to evaluate the magnitude of this eventual change, we shall use the same value aco = 0.65 mb at all 
energies. 
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Figure 1: Results on the centrality dependence of the J/ip nuclear modification in PbPb collisions 
at 2.76 TeV at mid (left) and forward (right) rapidity compared to ALICE data [351 EHl EZl EH]. 
The dashed line corresponds to the shadowing effect. The lowest line corresponds to the combined 
effect of shadowing + comover dissociation. The shadowed area corresponds to our result when 
shadowing + comover dissociation + recombination are taken into account. The uncertainty 
takes into account the variation between the minimum and maximal value of C. 



comovers, which induces a decrease of the C-value. This will produce a decrease of R/j^ with increasing 
y. Note nevertheless that this effect may be compensated by the increase of R^/^ due to a smaller 
density of comovers at y 7^ 0, which induces less dissociation. We have chosen a smooth behaviour of 
C{y) with rapidity, which in the rapidity range 2.5 < y < 4 corresponds to a mean reduction in the 
C-value of the order of 20%. Taking j*^ in this forward rapidity range, = 2.23 fxh ^2], these 

C-values correspond to an input charm cross section ~ 0.4-^0.6 mb in the forward rapidity region. 
Our procedure gives a reasonable description of data when recombination is taken into account. 
Some considerations proceed here: 

First, the behaviour of the nuclear modification factor i?^'^ for different rapidity ranges changes 
depending of the amount of recombination considered. If the recombination C-value is considered to 
be small -lower limit-, the total suppression will be controlled essentially by the comover dissociation, 
and the nuclear modification factor will follow the same behaviour as the usual comover suppression, 
but with a higher absolute value. On the other hand, when the C-value is taken as its upper limit. 
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the recombination term controls the total suppression, inducing a decrease of the nuclear modification 
factor when going to forward rapidities. 

Second, the shadowing corrections for heavy quarks we use here have been calculated in a very 
simple way, according to eq. Q within the Glauber-Gribov theory. The advantage is that all calcu- 
lations can be easily done analytically. This shadowing roughly agrees with EKS98/nDSg predictions 
\39\ \^3\ H5] , in particular in the mid-rapidity region. In Fig. [2] we show the comparison of 
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Figure 2: Centrality dependence of the Glauber-Gribov shadowing (eq. (4)) compared to EKS98/nDSg 
calculations performed according to |39] in PbPb collisions at ^/s = 2.76 TeV. 



the shadowing model applied here with some recent calculations developed within the framework of 
Refs. |44| I39j . One can see that the Glauber-Gribov inspired model agrees with the lowest band of 
nDSg according to [lH |39]. On the other side, calculations of the shadowing developed within the 
framework of the Color Evaporation model j45j , lead in general to an slightly higher suppression than 
|44j . Because of this we consider that the shadowing presented here agrees with the models men- 
tioned above within uncertainties. Nevertheless, the Glauber-Gribov shadowing is almost constant 
with rapidity. Even if the y-dependence of the different shadowing models in AA collisions and in the 
rapidity ranges here considered is quite flat, as it shown in [35] and Refs. \39\ I45j therein, we are aware 
that our approach can induce an small overestimation of the shadowing suppression in the forward 
rapidity region. In fact we contemplate the upgrading of the CIM by the introduction of EKS98/nDSg 
according to |39l Hi! ■ 
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We proceed now to study the behaviour of the nuclear modification factor when different transverse 
momentum cuts are introduced. In Fig. [s] (left) our results for 2 different px ranges, < < 2 GeV 
and 5 < < 8 GeV, are compared to ALICE forward rapidity data [38j. Clearly, the amount of 




Figure 3: Results on the centrality dependence of the J nuclear modification factor in PbPb colli- 
sions at 2.76 TeV compared to ALICE forward rapidity data [38] in the pT ranges < < 2 GeV 
and 5 < < 8 GeV (left) and compared to CMS mid-rapidity data [171 SSI US] in the pT range 
6.5 < < 30 GeV (right). The dashed line corresponds to the shadowing effect. The lowest line 
(right) corresponds to the combined effect of shadowing -|- comover dissociation. The shadowed 
area corresponds to our result when shadowing -|- comover dissociation -|- recombination are 
taken into account. The uncertainty takes into account the variation between the minimum and 
maximal value of C. 

recombination is more important in the low pT region. The comparison to CMS data at mid-rapidity 
and higher px, 6.5 < < 30 GeV [13 HHl HSj, in Fig. [s] (right) emphasizes this finding: the amount 
of recombination is much smaller than in the mid-rapidity range at low pT, as shown in Fig. [T] (left). 
It is important to remark that this experimental fact -stronger suppression at higher pT~ can not be 
due to initial-state shadowing effects: the presence of shadowing corrections, more relevant at lower 
PTi conspirates in the opposite direction. 
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Moreover, our results on the dependence of the J /if) nuclear modification factor on the transverse 
momentum are shown in Fie;. [4] compared to ALICE [S^ and CMS [171 08] data. An overall agreement 
is obtained. 




p^. (GeV/c) p^. (GeV/c) 

Figure 4: Results on the transverse dependence of the J /tjj nuclear modification in PbPb collisions at 
2.76 TeV compared to CMS data [HlllH] at mid-rapidity (left) and to ALICE [38] and CMS data [57] 
at forward rapidity (left). The CMS data on the left part of the figure corresponds to the rapidity 
range 1.6 < \y\ < 2.4, while the ALICE data lie in the range 2.5 < y < 4. The dashed line corresponds 
to the shadowing effect. The lowest line corresponds to the combined effect of shadowing + 
comover dissociation. The shadowed area corresponds to our result when shadowing + comover 
dissociation + recombination are taken into account. The uncertainty takes into account the 
variation between the minimum and maximal value of C. 

We continue by showing our results versus rapidity. In Fig. [s] (left) we compare our results with 
CMS data [471 HH] in the mid-rapidity region. An overall agreement is obtained. The fact that we do 
not reproduce the detailed behaviour of ALICE data |37j . as can be seen in Fig. [s] (right), is extremely 
instructive. As we have mentioned above, we have chosen a very conservative behaviour of C{y) with 
rapidity, which leads to a mean reduction in the C-value of the order of 20% in the rapidity range 
2.5 < y < 4. The data in this region agrees with this choice in the interval 2.5 < y < 3.5, while for the 
most forward points an scenario with smaller amount of recombination or only dissociation without 
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recombination is clearly favoured. 




Figure 5: Results on the rapidity dependence of the J nuclear modification in PbPb collisions at 
2.76 TeV compared to CMS data [171 HE] at mid-rapidity (left) and to ALICE [37j at mid and forward 
rapidity (left). The lowest line corresponds to the combined effect of shadowing + comover disso- 
ciation. The shadowed area corresponds to our result when shadowing + comover dissociation 
+ recombination are taken into account. The uncertainty takes into account the variation between 
the minimum and maximal value of C. 



We finish by presenting our predictions for the LHC energy of 5.5 TeV. Here we have let vary our 
C-value between 2 and 4. Taking BRu x \y=Q w 350 nb [33] and app = 60 mb, these C-values 
correspond to an input charm cross section ~ 0.8 -^1.15 mb in the mid-rapidity region. In the 
rapidity range 2.5 < y < 4 we take a mean reduction in the C-value of the order of 20%, which 
corresponds to an input charm cross section ~ 0.6 ~ 0.9 mb when BRu x 250 nb |33| . 

Our results in the mid and forward rapidity ranges are plotted in Fig. [6] 

Note that in the previous work |16j stronger shadowing corrections -of around 20% larger than the 
present ones- for the heavy quark production were considered, which implied both less recombination 
effects and stronger total suppression. This demonstrates that the implementation of the initial-state 
effects is more relevant than what is currently assumed in the recombination approaches, since it 
affects the probability of regeneration. 
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Figure 6: Results on the centrality dependence of the J/ip nuclear modification in PbPb collisions 
at 5.5 TeV at mid (left) and forward (right) rapidity. The dashed line corresponds to the shad- 
owing effect. The lowest line corresponds to the combined effect of shadowing + comover dis- 
sociation. The shadowed area corresponds to our result when shadowing + comover dissoci- 
ation + recombination are taken into account. The uncertainty takes into account the varia- 
tion between the minimum and maximal value of C. For the mid-rapidity region, we have taken, 
from down to up, = 0.8,0.85,0.9,0.95,1.0,1.1,1.15 mb and for the forward rapidity region 

^ = 0.65,0.67,0.71,0.74,0.82,0.84,0.87,0.89 mb. 



4 Conclusion 

In this work we have studied the combined effect of J /ip dissociation and recombination of cc pairs 
into J /ip in the comover interaction model. This model does not assume thermal equilibrium of the 
matter produced in the collision and includes a comprehensive treatment of initial-state effects, such 
as shadowing. We estimate the magnitude of the recombination term from J /ip and open charm yields 
in pp collisions at LHC. Without any adjustable parameters, the centrality, transverse momentum and 
rapidity dependence of experimental data is reproduced. 

In our approach, the magnitude of the recombination effect is controlled by the total charm cross 
section in pp collisions. Note that, contrary to the results in [16J, where the combined effect of initial- 
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state shadowing and comover dissociation appeared to overcome the effect of parton recombination at 
5.5 TeV, we find here that the recombination effects are of crucial importance in PbPb colhsions at 
LHC energies and can dominate over the suppression, in agreement with |4 H 1401 H6] . The reason for 
this dissent is the fact that a different approach for the shadowing factor for heavy quark production 
was used in [16J, which minimaUzed the amount of recombination and led to an overestimation of the 
total suppression. 

Also, as a general remark, we are aware that the comover interaction model at these energies 
should not be considered to describe a final-state interaction at the hadronic level. Indeed, at small 
values of the proper time these comovers should be considered as a dense partonic medium. A large 
contribution to the comover interaction comes from the few first fm/c, where the system is in partonic 
or pre-hadronic stage. 
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